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Abstract—VLBI measurements of the absolute proper motions of 23 radio stars 
have been collected from published data. These are stars with maser emission, or very 
young stars, or asymptotic-giant-branch stars. By comparing these measurements with 
the stellar proper motions from the optical catalogs of the Hipparcos Celestial Reference 
Frame (HCRF), we have found the components of the residual rotation vector of this 
frame relative to the inertial coordinate system; {ujx,uJy,uJz) = (—0.39,—0.51,—1.25) ± 
(0.58,0.57,0.56) mas yr“^. Based on all the available data, we have determined new 
values of the components of the residual rotation vector for the optical realization of the 
HCRF relative to the inertial coordinate system; {uj,,.,ujy,ujz) = (—0.15,-1-0.24,—0.53) ± 
(0.11,0.10,0.13) mas yr“^. 

INTRODUCTION 

The present-day standard astronomical coordinate system, the International Celestial 
Reference System (ICRS), has been an official one since January 1, 1998, by the lAU 
decision. The directions of the principal plane and the coordinate origin at the epoch 
J2000.0 were hxed in the lAU decision. 

This system is realized as a catalog of the positions of 212 compact extragalactic radio 
sources uniformly distributed over the entire sky. The catalog was compiled from very 
long baseline interferometry (VLBI) observations, which is reflected in Ma et ah (1998). 
This specihc realization of the International Celestial Reference System is designated as 
ICRF (International Celestial Reference Frame). 

The current version of the catalog, ICRF2, already uses 295 most reliable radio sources, 
and the total number of objects in it is 3414 (Ma et ah 2009). The parameters of the third 
version, ICRF3, are now being discussed (Jacobs et al. 2014). The selection of candidates 
for ICRF3 is planned to be hnished by 2018 in order that the catalog be ready before the 
completion (by about 2021) of the GAIA optical space experiment. 

The hrst realization of the ICRS in the optical range was the Hipparcos (1997) 
Catalogue, which is usually designated as HCRF (Hipparcos Celestial Reference Frame). 
As was shown by Kovalevsky et al. (1997), the Hipparcos system was tied to the 
extragalactic reference frame with an accuracy of ±0.6 mas for the coordinates at the 
epoch 1991.25 and with an error of ±0.25 mas yr“^ for the rotation around three axes. 
The difficulty of tying Hipparcos to the ICRF lies in the fact that there are virtually no 


1 




common objects, because in the optical range the quasars turned out to be too faint to 
be directly observed in the Hipparcos experiment. Therefore, the main observing program 
for tying Hipparcos to the ICRF was the program of VLBI observations of 12 radio stars 
located near quasars (Lestrade 1999). The programs to determine the absolute proper 
motions of stars relative to galaxies and quasars (Kovalevsky et ah 1997) also served to 
check the tying. 

The application of various analysis techniques shows that there is a slight residual 
rotation of the HCRF relative to the inertial coordinate system with ujz ~ —0.4 mas yr“^ 
(Bobylev 2004a, 2004b). To check whether the HCRF is inertial, VLBI observations of 46 
radio stars are being performed (Boboltz et ah 2003, 2007) in the USA using the “Very 
Large Array” (VLA) in combination with the Pie Town antenna (New Mexico) at 8.4 
GHz. 

In this paper, we want to focus our attention on the VLBI observations of masers and 
Mira variables that are performed to study the Galaxy (Reid et ah 2014; Nakagawa et 
al. 2014). The VLBI observations of such sources carried out in the last 5-6 years show 
that the absolute proper motions of radio stars have a very high accuracy, 0.5-0.05 mas 
yr“^. This accuracy has been achieved, first, owing to the long interferometer baselines 
and, second, owing to the observations at high frequencies, 22.2 GHz and even 43.2 GHz 
(SiO masers). The goal of this work is to study the possibility of using such observations 
to determine the residual rotation of the HGRF relative to the present-day realization of 
the inertial coordinate system. 

VLBI DATA 

One of the projects to measure the trigonometric parallaxes and proper motions is 
the Japanese VERA (VLBI Exploration of Radio Astrometry) project devoted to the 
observations of H 2 O masers at 22.2 GHz (Hirota et ah 2007) and SiO masers at 43.2 GHz 
(Kim et al. 2008). 

Methanol (GH 3 OH, 6.7 and 12.2 GHz) and H 2 O masers are being observed in the USA 
on the VLB A (Reid et al. 2014). Similar observations are also being carried out within the 
framework of the European VLBI network (Rygl et al. 2010). These two programs enter 
into the general BeSSeL0 (Bar and Spiral Structure Legacy Survey) project (Brunthaler 
et al. 2011 ). 

For the same purpose, the VLBI observations of radio stars are being carried out in 
continuum at 8.4 GHz (Torres et al. 2007; Dzib et al. 2012). The radio sources located in 
the local arm (Orion arm) that are associated with young low-mass protostars are being 
observed within the framework of this program. 

Mira variables at the asymptotic-giant-branch (AGB) stage are being observed in the 
VERA program. Such stars manifest themselves in the radio band as masers (Nakagawa et 
al. 2014). Thus, our sample includes young massive supergiants, intermediate-mass giants, 
and low-mass T Tauri stars. 

An important peculiarity of all the listed observations is that the stellar positions have 
always been determined using two or three extragalactic sources, i.e., the parallaxes and 

^http://www3.mpifr-bonn.mpg.de/staff/abrunthaler/BeSSeL/index.shtml 
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TaGjilma 1: Characteristics of the radio stars 


Star 

HIP/UCAC 

Type 

Erequency, 

GHz 

Program 

Emission 

Ref 

T Lep 

HIP 23636 

Mira 

22.2 

VERA 

H 2 O masers 

(1) 

S Crt 

HIP 57917 

SRb 

22.2 

VERA 

H 2 O masers 

(2) 

W Hya 

HIP 67419 

SRa 

1.6 

NRAO 

OH masers 

(3) 

RX Boo 

HIP 70401 

SRb 

22.2 

VERA 

H 2 O masers 

(4) 

S CrB 

HIP 75143 

Mira 

1.6 

NRAO 

OH masers 

(5) 

U Her 

HIP 80488 

Mira 

1.6 

NRAO 

OH masers 

(5) 

RR Aql 

HIP 98220 

Mira 

1.6 

NRAO 

OH masers 

(5) 

R Aqr 

HIP117054 

Mira 

43.2 

VERA 

SiO masers 

(6) 

R Cas 

HIP118188 

Mira 

1.6 

NRAO 

OH masers 

(3) 

SY Scl 

UCAC4 

Mira 

22.2 

VERA 

H 2 O masers 

(7) 

UX Cyg 

UCAC4 

Mira 

22.2 

NRAO 

H 2 O masers 

(8) 

ss Cyg 

UCAC4 

D. Nova 

8.4 

NRAO 

continuum 

(9) 

IM Peg 

HIP112997 

RS CVn 

8.4 

VLBI 

continuum 

(10) 

S Per 

HIP 11093 

SRc 

22.2 

NRAO 

H 2 O masers 

(11) 

V773 Tau 

HIP 19762 

T Tau 

8.4 

NRAO 

continuum 

(12) 

HDE 283572 

HIP 20388 

T Tau 

8.4 

NRAO 

continuum 

(13) 

T Tau N 

HIP 20390 

T Tau 

8.4 

NRAO 

continuum 

(14) 

LSI +61 303 

HIP 12469 

XMXRB 

8.4 

VLBI 

continuum 

(15) 

Cyg X-1 

HIP 98298 

XMXRB 

8.4 

NRAO 

continuum 

(16) 

Cyg OB2#5 

HIP101341 

EB 

8.4 

NRAO 

continuum 

(17) 

IRAS 22480+6002 

UCAC4 

— 

22.2 

VERA 

H 2 O masers 

(18) 

PZ Cas 

HIP 117078 

SRa 

22.2 

VERA 

H 2 O masers 

(19) 

9 ^ Ori A 

UCAC4 

— 

8.4 

NRAO 

continuum 

(20) 


SR—semiregular pulsating stars; RS—eruptive variables of the RS Canum Venaticorum type; 
EB—eclipsing binaries; XMXRB—high-mass X-ray binaries; D. Nova—dwarf novae. 

(1): Nakagawa et al. (2014), (2): Nakagawa et al. (2008), (3): Vlemmings et al. (2003), (4): 
Kamezaki et al. (2012), (5): Vlemmings and van Langevelde (2007), (6): Min et al. (2014), (7): 
Nyu et al. (2011), (8): Kurayama et al. (2005), (9): Miller-Jones et al. (2013), (10): Ratner et 
al. (2012), (11): Asaki et al. (2010), (12): Torres et al. (2012), (13): Torres et al. (2007), (14): 
Loinard et al. (2007), (15): Dhawan et al. (2006), (16): Reid et al. (2011), (17): Dzib et al. (2013), 
(18): Imai et al. (2012), (19): Kusuno et al. (2013), (20): Menten et al. (2007). 
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Ta6jiHii,a 2; Proper motions of the radio stars 



VLBI 

HIPPARCOS-2007 

Star 

Ha COS (5 ± e, 
mas yr“^ 

Hs + e, 
mas yr“^ 

Ha COS 5 ±e, 
mas yr“^ 

Hs + e, 
mas yr“^ 

T Lep 

14.60 + 0.50 

-35.43 + 0.70 

11.43 + 0.98 

-33.34 + 0.86 

S Crt 

-3.17 + 0.22 

-5.41 + 0.22 

-3.37 + 1.00 

-4.68 + 0.75 

W Hya 

-44.24 + 2.04 

-55.28 + 2.98 

-49.31 + 1.48 

-59.71 + 1.04 

RX Boo 

24.55 + 1.06 

-49.67 + 2.38 

21.21 + 0.50 

-48.79 + 0.46 

S CrB 

-9.06 + 0.23 

-12.52 + 0.29 

-7.73 + 0.57 

-13.03 + 1.02 

U Her 

-14.98 + 0.29 

-9.23 + 0.32 

-15.94 + 0.61 

-11.03 + 0.69 

RR Aql 

-25.11 + 0.74 

-49.82 + 0.54 

-25.12 + 4.46 

-49.14 + 3.23 

R Aqr 

37.13 + 0.47 

-28.62 + 0.44 

33.00 + 1.53 

-25.74 + 1.30 

R Cas 

80.52 + 2.35 

17.10 + 1.75 

85.52 + 0.75 

17.49 + 0.72 

SY Scl 

5.57 + 0.04 

-7.32 + 0.12 

4.60 + 0.90 

-7.80 + 0.90 

UX Cyg 

-6.91 + 0.75 

-12.52 + 1.57 

-8.90 + 1.60 

-9.80 + 2.50 

ss Cyg 

112.42 + 0.07 

33.38 + 0.07 

113.20 + 0.90 

33.40 + 1.00 

IM Peg 

-20.83 + 0.09 

-27.27 + 0.09 

-20.73 + 0.28 

-27.75 + 0.27 

S Per 

-0.49 + 0.35 

-1.19 + 0.33 

-2.70 + 2.20 

-0.29 + 1.65 

V773 Tau 

8.30 + 0.50 

-23.60 + 0.50 

4.11 + 2.74 

-24.48 + 1.88 

HDE 283572 

8.88 + 0.06 

-26.60 + 0.10 

6.84 + 1.64 

-27.15 + 1.12 

T Tau N 

12.35 + 0.04 

-12.80 + 0.05 

15.51 + 1.93 

-13.67 + 1.64 

LSI +61 303 

-0.30 + 0.07 

-0.26 + 0.05 

0.27 + 2.91 

2.38 + 2.52 

Cyg X-1 

-3.78 + 0.06 

-6.40 + 0.12 

-3.37 + 0.75 

-7.15 + 0.86 

Cyg OB2#5 

-1.64 + 0.98 

-7.16 + 1.28 

-2.14 + 1.23 

-2.48 + 1.14 

IRAS 22480+6002 

-2.58 + 0.33 

-1.91 + 0.17 

-2.60 + 1.90 

-2.10 + 6.10 

PZ Cas 

-3.70 + 0.20 

-2.00 + 0.30 

-4.15 + 0.85 

-3.55 + 0.81 

6^ Ori A 

4.82 + 0.09 

-1.54 + 0.18 

-0.90 + 2.50 

0.10 + 2.50 


proper motions of these radio stars are absolute. 

Table 1 presents some characteristics of the radio stars. Columns 1, 2, 3, 4, 5, 6, 
and 7 give, respectively, the star name, the number in the Hipparcos Catalogue, the 
type of variability, the observation frequency, the name of the observing program or 
observatory, the type of emission characterizing the VLBI observations (for example, they 
were performed in maser lines or in continuum), and references. Note that the Hipparcos 
proper motions are very unreliable for the massive O star 6^ Ori A (HIP 26220), a member 
of the famous Orion Trapezium; therefore, we used the data from the UCAC4 (2012) 
catalog. It can be seen from the second column of Table 1 that there are four more stars 
(absent in the Hipparcos Catalogue) with the UCAC4 proper motions. 

We did not include the red supergiant VY CMa (HIP 35793) in our list, for which high- 
accuracy VLBI observations are available (Zhang et ah 2012). For this star, the absolute 
values of the “Hipparcos minus VLBI” proper motion differences are about 8 mas yr“^ in 
each of the coordinates. Such a large difference stems from the fact that slightly different 
parts of the extended asymmetric envelope of this star are observed in the optical and 
radio bands. 
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Note that long-term observations in maser lines are difficult to carry out, because 
the maser spots are highly unstable; therefore, the periods of observations for such stars 
are 1.5-2 years (the minimum period required to determine the annual parallax). The 
situation is different for continuum observations. For example, the period of observations 
for IMPeg was about 8 years (Ratner et ah 2012). 

Table 2 presents the proper motions of the radio stars measured by the VLBI technique 
(according to the references in Table 1) and the proper motions of these stars from the 
version of the Hipparcos Catalogue revised by van Leeuwen (2007). Van Leeuwen (2007) 
showed that the new version of the Catalogue completely reproduces the previous system, 
i.e., it has no residual rotation relative to the Hipparcos-1997 version, while it excels 
considerably the Hipparcos-1997 version in terms of random errors (especially in the 
region of bright stars). It is of interest to check this using the available material. 

Figure 1 presents the “HIP-1997 minus VLBI” and “HIP-2007 minus VLBI” stellar 
proper motion differences. It can be clearly seen from our comparison of the graphs that 
the dispersion of the “HIP-1997 minus VLBI” differences exceeds considerably that of 
the “HIP-2007 minus VLBI” differences. For some of the stars, the absolute values of the 
“HIP-1997 minus VLBI” differences in right ascension exceed 6 mas yr“^. 

In Fig. 2, the “HIP-2007 minus VLBI” stellar proper motions are plotted against the 
equatorial coordinates a and S. It can be seen from the hgure that a wave in a and a 
slight trend in S are noticeable in the differences A/i^ cos <5, while the differences A/i^ are 
distributed quite symmetrically relative to the horizontal axis. 

The following coupling equations can be used to determine the three angular velocities 
of mutual rotation of the two systems Ux,ujy,Uz : 

Afia cos S = —Ux cos a sin 5 — Uy sin a sin 5 -|- Uz cos 5, . ^ 

A/i 5 = cja; sin a — cjy cosa, ^ 

where the “Hipparcos minus VLBI” differences are on the left-hand sides of the equations. 
Having solved the system of conditional equations (1) by the least squares method for 
the “HIP-1997 minus VLBI” difference, we obtained the following rotation components 
(in mas yr“^): 

Ux = —0.85 ± 0.63, 

=+0.05 ±0.61, (2) 

= -1.33 ±0.61, 

where the error per unit weight, which characterizes the dispersion of the residuals, is 
(To = 2.37 mas yr“^. For the “HIP-2007 minus VLBI” difference, we obtained the following 
rotation components (in mas yr“^); 

ujx = —0.39 ± 0.58, 

Wy =-0.51 ± 0.57, (3) 

= -1.25 ±0.56, 

where the error per unit weight is (Tq = 2.19 mas yr“^. We see that it is more advantageous 
to use the HIP-2007 version. 
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A|J„* masyr-l 

Phc. 1; (a) “HIP-1997 minus VLBI” stellar proper motion differences and (b) “HIP-2007 
minus VLBI” stellar proper motion differences; A/Tq* = Afi^cosS. 


HCRF INERTIALITY 

Table 3 presents almost all of the results of comparing the individual programs with the 
HCRF catalogs known to date. The first seven rows of Table 3 give the results used 
by Kovalevsky et al. (1997) to calibrate Hipparcos and to estimate its residual rotation 
relative to the system of extragalactic sources. 

(1) The “VLBI-1999” solution was obtained by comparing the absolute proper motions 
of twelve radio stars with the Hipparcos Catalogue. A detailed description of the 
observations can be found in Lestrade et al. (1999). In total, 21 antennas of various 
diameters in the USA and Europe were involved in the observations. The period of 
observations for each star was from 2 to 11 years. The observations were performed at 5.0 
and 8.4 GHz. The coordinates, parallaxes, proper motions, and even accelerations in the 
proper motions of stars were determined. According to the estimates by Lestrade et al. 
(1999), the error in the stellar proper motion for some of the stars was about 0.05 mas 
yr-P 

(2) As the “NPMl” solution, we used the results of comparing the stellar proper 
motions from the NPMl (Klemola et al. 1994) and Hipparcos Catalogues by the Heidelberg 
team. We took the parameters of this solution from Table 2 in Kovalevsky et al. (1997). 
This solution was obtained in the range of magnitudes 10"*.5-11™'.5, where (Fig. 1 in 
Platais et al. (1998a)) the “HIP minus NPMl” stellar proper motion differences have a 
“horizontal” character, near zero. In our opinion, the proper motions from the NPMl 
Catalogue in this magnitude range are free to the greatest extent from the magnitude 
equation, which is considerable in this Catalogue. 

(3) The “KIEV” solution. The stellar proper motions from theGPMl(Rybka and 
Yatsenko 1997) and Hipparcos Gatalogues were compared by Kislyuk et al. (1997). 
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Phc. 2: “HIP-1997 minus VLBI” stellar proper motion differences versus equatorial 
coordinates a and 5; A/Tq,* = AfiaCosS. 


(4) The “POTSDAM” solution was taken from Hirte et al. (1996). 

(5) The “BONN” solution was obtained by performing the Bonn program, which is 
reflected in Geffert et al. (1997) and Tucholke et al. (1997). 

(6) The “EOP” solution was obtained as a result of the analysis of the Earth Orientation 
Parameters (EOPs) by Vondrak et al. (1997). Only two orientation parameters, and 
Uy, are determined in this method. 

(7) The “HST” solution was obtained by Hemenway et al. (1997) using the Hubble 
Space Telescope (HST). Note that this solution has virtually no influence on the 
calculation of the weighted average because of the large random errors in the comparison 
parameters. 

(8) The “SPM2” solution was obtained by Zhu (2001) by comparing the stellar proper 
motions from the SPM2 (Platais et al. 1998b) and Hipparcos Catalogues. 

(9) The “PUL2” solution was found by comparing the PUL2 Pulkovo photographic 
catalog (Bobylev et al. 2004) and the Hipparcos Catalogue. 

(10) The “XPM” solution is based on the Khar’kov catalog of absolute stellar proper 
motions, XPM (Fedorov et al. 2009). It was absolutized using ?^1.5 million galaxies from 
the 2MASS catalog of extended sources (Skrutskie et al. 2006). Thus, the XPM catalog 
is an independent realization of the inertial coordinate system. The proper motions from 
the XPM and UCAC2 (Zacharias et al. 2004) catalogs were compared in Bobylev et al. 
(2010), where the parameters Ux, ooy, and Uz were calculated using f»l million stars. 

(11) The “MINOR PLANETS” solution. Chernetenko et al. (2008) estimated the 
orientation parameters of the Hipparcos system relative to the coordinate systems of the 
DE403 and DE405 ephemerides by analyzing a long-term series of asteroid observations. 
This result leads to the conclusion that either the DE403 and DE405 dynamical theories 
need an improvement or the Hipparcos system needs a correction. We reduced the weight 
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of this solution by half because of the possible contribution from the inaccuracy of the 
DE403 and DE405 dynamical theories. 

(12) The “VLA + PT-2007” solution was obtained by Boboltz et al. (2007) by analyzing 
the positions and proper motions of 46 radio stars. The VLBI observations were performed 
in the USA using the “Very Large Array” (VLA) in combination with the Pie Town antenna 
(New Mexico) at 8.4 GHz. In comparison with the observations by Lestrade et al. (1999), 
these observations have a considerably lower resolution in positional observations, because 
the interferometer size is smaller. On average, the positional errors are 13 mas in right 
ascension and 16 mas in declination (Boboltz et al. 2003). To achieve a high accuracy 
of determining the proper motions, they used long series of observations, about 20 years 
for each star. Only the coordinates and proper motions of radio stars were the quantities 
being determined. According to the estimates by Boboltz et al. (2007), the error in the 
stellar proper motion in this program is, on average, f^l.7 mas yr“^ in each coordinate. 

(13) The search for the “VLBI-2014” solution is described in the first part of this paper. 
Note that the three solutions under consideration based on the VLBI observations of radio 
stars, “VLBI-1999”, “VLA + PT-2007”, and “VLBI-2014”, have significant differences in 
observing techniques; therefore, we consider them as independent solutions. 

The weight assigned to each comparison catalog is inversely proportional to the square 
of the mean error in the corresponding quantities (ux,u)y,Uz) and was calculated from 
the formula 

Pi = efc^e„VG^ i = 1, •••, 13. (4) 

Equations of the form (1), where the “Hipparcos minus Catalogue” differences are on the 
left-hand sides, can serve to determine Ux, cOy-, and Uz- Lindegren and Kovalevsky (1995) 
proposed a slightly different form of these equations; 

Afia COS 6 = Ux cos a sin (5 -|- Uy sin a sin S — Uz cos 6, , . 

A/i^ = —cUa; sin a 4 -cosa, ^ 

where the “Catalogue minus Hipparcos” differences are on the left-hand sides of the 
equations. Zhu (2001) and Boboltz et al. (2007) published the parameters that they 
determined with a change in the form of either the left-hand sides or the righthand sides 
of Eqs. (1) or (5). In these two cases, we brought the signs of the quoted parameters to 
the necessary standard form. 

The last rows of Table 3 present “Average 1” that was calculated as a simple average 
but without using the HST solution and “Average 2” that was calculated from all data 
as a weighted mean and is the main result of our analysis. It is important to note that 
the rotation component Uz = —0.53 ± 0.13 mas yr“^ in the solution obtained differs 
signihcantly from zero. 

DISCUSSION 

By comparing the XPM stars with catalogs extending the HCRF to faint magnitudes, 
such as PPMXL (Roeser et al. 2010), UCAC3 (Zacharias et al. 2009), Tycho-2 (Hog et al. 
2000), and XCl (Fedorov and Myznikov 2006), Fedorov et al. (2011) showed that the two 



Ta6jiima 3; Components of the residual rotation vector of the optical realization of the 
ICRS/Hipparcos system relative to the inertial coordinate system 


Method 

N. 

N 

’ area 

ujx, mas yr ^ 

Uy, mas yr ^ 

cuz, mas yr ^ 

VLBI-1999 

12 


-0.16 ±0.30 

-0.17 ±0.26 

-0.33 ±0.30 

NPMl 

2616 

899 

-0.76 ±0.25 

±0.17 ±0.20 

-0.85 ±0.20 

Kiev 

415 

154 

-0.27 ±0.80 

±0.15 ±0.60 

-1.07 ±0.80 

Potsdam 

256 

24 

±0.22 ±0.52 

±0.43 ±0.50 

±0.13 ±0.48 

Bonn 

88 

13 

±0.16 ±0.34 

-0.32 ±0.25 

±0.17 ±0.33 

EOP 



-0.93 ±0.28 

-0.32 ±0.28 

— 

HST 

78 


-1.60 ±2.87 

-1.92 ± 1.54 

±2.26 ±3.42 

SPM2 

9356 

156 

±0.10 ±0.17 

±0.48 ±0.14 

-0.17 ±0.15 

PUL2 

1004 

147 

-0.98 ±0.47 

-0.03 ±0.38 

-1.66 ±0.42 

XPM 

1 X 10® 

1431 

-0.06 ±0.15 

±0.17 ±0.14 

-0.84 ±0.14 

Minor Planets 

116 


±0.12 ±0.08 

±0.66 ±0.09 

-0.56 ±0.16 

VLA±PT-2007 

46 


-0.55 ±0.34 

-0.02 ±0.36 

±0.41 ±0.37 

VLBI-2014 

23 


-0.39 ±0.58 

-0.51 ±0.57 

-1.25 ±0.56 

Average 1 



-0.29 ±0.12 

±0.06 ±0.10 

-0.55 ±0.20 

Average 2 



-0.15 ±0.11 

±0.24 ±0.10 

-0.53 ±0.13 


is the number of stars/asteroids, is the number of areas on the celestial sphere, average 1 
is a simple average (without HST); average 2 is a weighted average. 


rotation components Ux and uiy do not differ significantly from zero, while the Z rotation 
component has a measurable value, Uz = —1.8 ± 0.16 mas yr“^. 

Note the result by Bobylev and Khovritchev (2011) obtained from the “proper motions” 
of about 8000 galaxies from the UCAC3 catalog (Zacharias et al. 2009). These authors 
found the residual rotation of the HCRF with respect to the inertial coordinate system 
(rij., r2;j) = (0.58,-1.02,-0.59) ± (0.15,0.15,0.17) mas yr“^ in Galactic coordinates. 
The corresponding rotation components around the equatorial axes (in mas yr“^) are 

Ux = -0.02 ±0.15, 

cuy = ±0.06 ±0.15, (6) 

= -1.31 ±0.17. 

Fedorov et al. (2011), Wu et al. (2011), and Grabowski et al. (2014) showed the “proper 
motions” of galaxies and quasars from the PPMXL catalogs to have measurable values. 
The proper motions averaged over all galaxies and quasars are ~ —2 mas yr“^ in each of 
the coordinates {fiaCosS and fis)- Wu et al. (2011) and Grabowski et al. (2014) did not 
determine the global rotation parameters, but they proposed to perform an additional 
absolute calibration of the stellar proper motions from the PPMXL catalog using a table 
of corrections. 

Another interesting example of allowance for the additional absolute calibration of 
the stellar proper motions from the PPMXL catalog using galaxies is given in Lopez- 
Gorredoira et al. (2014). These authors studied the parameters of such a large-scale 
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phenomenon as the Galactic warp using red giant clump stars. It can be clearly seen 
from Fig. 3 in the cited paper that the corrections to the vertical stellar velocities reach 
about 70 km s“^ at a heliocentric distance of about 8 kpc toward the Galactic anticenter. 
In other words, this distribution can be interpreted as the existence of residual rotation 
of the PPMXL stars around the Galactic Y axis with an angular velocity Qy ^ —1.5 mas 
yr“^, in good agreement with the result of analyzing the PPMXL catalog obtained by 
Fedorov et ah (2011). 

All of the aforesaid leads us to conclude that although the values of Gz that we derived 
from radio stars in solutions (2) and (3) have fairly large errors, they are in agreement with 
the results of the analysis of large catalogs with galaxies presented above, in particular, 
with the result ( 6 ) and the result from Fedorov et al. (2011). 

Having analyzed the individual solutions but using a smaller amount of data, Bobylev 
(2010) found the following rotation components: {ux,ujy,ujz) = (—0.11,+0.24,—0.52) ± 
(0.14, 0.10, 0.16) mas yr“^. The advantage of the result of this study obtained in a similar 
way (“Average 2” in Table 3) is that the random errors in the parameters decreased. Thus, 
each result of the individual comparison of catalogs is of great importance for the solution 
of our problem. 

CONCLUSIONS 

We collected highly accurate VLBI determinations of the absolute proper motions of 23 
radio stars from published data. The observations of the parallaxes and proper motions 
of these stars have been performed by various scientific teams in the last five or six years. 
The stars have a different evolutionary status. Some of them are very young stars with 
maser emission (H 2 O and GH 3 OH masers). AGB stars observed as OH, H 2 O, and SiO 
masers constitute the other part of the sample. Several stars were observed in continuum. 

By comparing these measurements with the proper motions from the HGRF catalogs, 
we found the components of the residual rotation vector of this system with respect 
to the present-day realization of the inertial coordinate system, namely: {ux,u}y,Uz) = 
(—0.39, —0.51, —1.25)±(0.58, 0.57, 0.56) mas yr“L This estimate is a completely new result 
of the individual comparison of the sample of radio stars with the Hipparcos Gatalogue. 

Based on all the available individual results, we determined new, most probable values 
of the components of the residual rotation vector for the optical realization of the HGRF 
with respect to the inertial coordinate system: {ux,ujy,ujz) = (—0.15,+0.24,—0.53) ± 
(0.11,0.10,0.13) mas yr“^. 
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